UNCLASSIFIED 
AD  NUMBER 


AD008032 

CLASSIFICATION  CHANGES 

TO: 

unclassified 

FROM: 

confidential 

LIMITATION  CHANGES 

TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Controlling  DoD  Organization:  Office  of  Naval 
Research,  Arlington,  VA  22203. 


AUTHORITY 

ONR  ltr  dtd  18  May  1966;  ONR  ltr  dtd  26  Oct 
1977 


THIS  PAGE  IS  UNCLASSIFIED 


UNCLASSIFIED 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 


CAMERON  STATION  ALEXANDRIA.  VIRGINIA 


Mi?i  i'.r 


RACED  AT  3 YEAR  INTERVALS: 
DECLASSIFIED  AFTER  12  YEARS 
DOD  DIR  5200  \0 


UNCLASSIFIED 


THIS  REPORT  HAS  BEEN  DECLASSIFIED 
AND  CLEARED  FOR  PUBLIC  RELEASE. 


DISTRIBUTION  A 
APPROVED  FOR  PUBLIC  RELEASE ; 
DISTRIBUTION  UNLIMITED. 


ONFIDENTIAL 


VELOCITY  OF  DETONATION  IN  'IDEAL' 
EXPLOSIVES  WITH  INERT  ADDITIVES 


TECHNICAL  REPORT 
NUMBER  X 
MARCH  25,1953 


CONTRACT  NUMBER 
N7  - onr  - 45 107 
PROJECT  NUMBER 
357  - 239 


EXPLOSIVES  RESEARCH  GROUP 
VSTITUTE  FOR  THE  STUDY  OF  RATE  PROCESSES 
UNIVERSITY  OF  UTAH 
SALT  LAKE  C I TY 


COPY  NO.  36 


CONFIDENTIAL 

SECURITY  IKFORHATI OK 


v 


Velocity  of  Detonation  in 
'Ideal'  Explosives  with  Inert  Additives 


Investigators:  Report  Prepared  by: 

M.  A.  Cock  D.  W.  Robinson 

D.  W.  Robinson 

Project  Director: 

M.  A.  Cook 

Technical  Report  No.  X 
March  25,  1953 


Contract  No.  N7-onr-45107 
Project  No.  357  239 


"THIS  MATERIAL  CONTAINS  INFORMATION  AFFECTING  THE  NATIONAL 
DEFENSE  OF  THE  UNITED  STATES  WITHIN  THE  MEANING  OF  THE 
ESPIONAGE  LAWS,  TITLE  18,  U.  S.  C. , SECTIONS  793  and  794 
THE  TRANSMISSION  OR  REVELATION  OF  WHICH  IN  ANY  MANNER  TO 
AN  UNUTHORIZED  PERSON  13  PROHIBITED  BY  LAW." 


\ 

Explosives  Research  Group 
Institute  for  the  Study  of  Rate  Processes 
University  of  Utah 
Salt  Lake  City 


CONFIDENTIAL 
SECURITY  INFORMATION 


i 


CONFIDENTIAL 


Velocity  of  Detonation  in 
•Ideal'  Explosives  with  Inert  Additivco 


Abstract 

A simplified  approximate  method  for  estimating  the  'ideal*  detonation 
velocity  of  pure  explosives  is  given.  Theoretical  applications  of  this 
method  are  included  to  study  the  effect  upcn  velocity  of  initial  tempera- 
ture and  of  incomplete  reaction  of  the  explosive.  A method  for  calculating 
the  ideal  velocity  cf  a pure  explosive  with  inert  ingredients  is  also  given. 
Several  different  mixtures  a^c  considered,  and  the  calculated  velocities 
are  reported.  Where  experimental  data  are  available,  the  calculated  and 
observed  results  are  compared. 

Introduction 

Solution  of  the  thermo-hydrodynamic  problem  for  an  explosive  by  the 
"inverse"  method  is  dependent  upon  an  experimental  'ideal'  detonation 
velocity  vs.  initial  density  relationship.  (An  'ideal'  detonating  explo- 
sive is  defined  as  one  in  which  the  velocity  does  not  increase  with  in- 
creasing diameter  of  the  charge.)  Using  this  empirical  relationship,  the 
various  thermodynamic  variables  for  the  detonation  proce  s can  be  cal- 
culated. Difficulty  immediately  arises,  however,  if  such  experimental 
data  for  a particular  explosive  are  not  available.  No  completely  satis- 
factory method  for  predicting  detonation  velocities  has  yet  been  suggested. 
The  first  part  of  this  report  is  ar.  attempt  to  present  a simplified  method 
for  predicting  approximately  the  ideal  velocities  of  pure  explosives,  hence, 
providing  an  approximate  solution  to  the  thermo-hydrodynamic  problem.  Due 
to  the  assumptions  which  are  made  in  this  method,  only  an  estimate  of  the 
velocities  can  be  expected.  The  method  was  applied,  however,  to  some 
explosives  with  known  velocities,  and  the  results  agreed  favorably  with 
the  observed  values. 

An  investigation  was  also  made  of  the  effect  of  solid  diluents  in 
pu^e  explosives.  The  second  part  of  this  paper  developes  a method  by 
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which  velocities  may  be  predicted  approximately  for  explosive-inert  mix- 
tures. This  method  was  successfully  applied  in  calculations  of  the 
detonation  velocities  for  several  mixtures,  agreement  v 1th  the  observed 
values  being  usually  within  experimental  error  of  measurement.  This 
theory  includes  the  compressibility  of  the  inert  substance.  Hence,  if 
measured  velocities  of  a particular  mixture  were  given,  by  an  '‘inverse" 
solution  the  compressibility  of  the  inert  could  be  calculated. 

Ideal  Velocities  in  Pure  Explosives 

It  was  demonstrated  in  T.  R.  VI  that  any  reasonable  choice  of  the 
form  of  the  equation  of  state  may  be  used  with  satisfactory  results  in 
the  solution  of  the  thermo-hydrodynamic  problem  of  an  explosive.  In 
this  report  the  relation 

pv  = nRTeX  (l) 

or  equivalently,  the  equation 

pv  = nRT  + pa  (2) 

was  used  as  the  form  of  the  equation  of  state.  By  eliminating  between 
(l)  and  (2),  the  relation 

a » v(l  - e-x)  (3) 

is  obtained.  Equation  (3)  explicitly  relates  the  two  forms  (l)  and  (2). 
Furthermore,  it  was  concluded  in  T.  R.  VI  that  present  methods  of  measur- 
ing thermo-hyarodynamic  variables  are  insufficient  to  distinguish  between 
reasonable  forms  of  equations  of  state.  The  following  is  an  attempt  to 
consider  the  converse  of  this  conclusion.  That  is,  by  "fixing"  an  equa- 
tion of  state,  assumed  valid  for  most  explosives,  it  should  be  possible 
to  make  reasonable  theoretical  estimates  for  the  values  of  the  thermo- 
hydrodynamic quantities,  including  detonation  velocity.  Since  the  "in- 
verse" method  of  solving  the  thermo-hydrodynamic  problem  requires  only 
one  empirical  relationship,  namely  the  velocity-density  curve,  it  is 
necessary  in  the  proposed  method  orfy  to  reduce  the  number  of  degrees  of 

• >rf*rPT  a t 
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freedom  of  (l)  or  (2)  by  one,  in  order  that  the  proposal  be  realized. 

The  fundamental  assumption  made  to  satisfy  this  condition  is  that  the 
variable  x and  equivalently  a are  known  functions  of  v.  That  is,  there 
exist  'universal'  relationships  x --  X(  v^  ) and  a = A(v2),  valid  for 
(most)  explosives.  This  assumption  is  not  justified  explicitly,  but  the 
following  comments  are  suggested  in  respect  to  the  particular  choice: 

(1)  There  is  a need  for  a simple,  practical  methcd  of  predicting 
velocities.  The  above  fundamental  assumption  leads  to  such  a method. 

The  results,  dependent  on  the  validity  of  the  assumption,  can  easily  be 
checked  with  some  known  calculations  or  observed  data. 

(2)  The  fundamental  assumption  agrees  with  a previous  observation 
of  Cook^;  viz,  that  most  explosive^  satisfy  the  same  a vs.  V2  relation- 
ship. Figure  1 is  a reprint  of  the  results  of  his  study. 

(3)  In  T.  k.  VI  two  different  sets  of  experimental  velocity  data 
were  used  in  determining  an  equation  of  state  for  EDNA.  If  it  were 
assumed  that  these  two  sets  of  data  represented  the  extreme  possibilities, 
then  the  effect  of  experimental  error  of  velocity  measurement  or.  the 
equation  of  state  would  be  that  indicated  in  figure  2.  Here  a plot  of 

x vs.  V2  is  shown  obtained  from  the  results  of  the  two  sets  of  calcula- 
tions given  in  tables  IV  and  V of  T.  R.  VI.  It  is  clear  that  any  choice 
of  a functional  relationship  x = X(v^) , which  is  within  the  limits  given 
here,  will  predict  the  velocity  of  EDNA  within  the  velocity  limits  of  the 
two  sets  of  results.  It  is  also  to  be  noted  that  there  exists  a liberal 
choice  for  such  a function  for  this  particular  explosive.  And  under  the 
assumptions  stated,  it  is  possible  to  accept  the  results  of  calculations 
obtained  by  use  of  any  such  function.  Thus,  it  is  reasonable  to  assume 
from  these  results  that  the  choice  of  a universal  relationship  x * 
will  not  be  too  critical  if  one  is  satisfied  to  compute  velocities  which 
are  accurate  only  to  within  about  + 150  m/sec. 

(U)  After  a specific  choice,  satisfying  the  fundamental  assumption, 
has  been  made,  the  selection  may  or  may  not  be  satisfactory.  If  the 
latter  is  the  case,  then  one  can  assume  that  an  incorrect  selection  has 
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been  na.de,  and  for  practical  purposes  attempt  to  correct  or  try  to  classify 
the  explosives  such  that  certain  relationships  will  be  valid  for  particu- 
lar classes.  Such  work  will  depend  upon  further  experimental  investiga- 
tions and  theoretical  calculations. 

If  the  scheme  is  successful,  it  provides  a simple  method  for  deter- 
mining velocities.  The  reasons  for  accepting  the  above  fundamental 
assumption  can  then  be  sought  with  profit. 

For  an  actual  selection  of  a functional  relationship  x * X(v2),  for 
this  report  all  of  the  available  evidence  at  hand  was  used.  It  was  se- 
lected as  a near  approximation  to  the  data  of  T.  R.  VI,  VIII,  and 
"Theory  of  Detonation" . ^ The  plot  is  given  in  figure  3,  and  a set  of 
corresponding  values  is  given  in  table  I. 

Two  methods  of  solving  the  thermo-hydrodynamic  problem  will  now  be 
. described.  The  first  introduces  the  fundamental  assumption  which  has 

just  been  discussed.  The  second  considers  further  assumptions,  which 
v permits  a simpler  solution.  The  first  requires  a calculation  of  the 

products  of  detonation  over  the  complete  range  of  temperature  and  fuga- 
city  factors  of  interest  for  a particular  explosive  compound.  The 
method  given  in  the  first  part  of  T.  R.  VIII  i3  available  for  this  cal- 
culation. With  this  known  composition,  the  total  number  of  n of  moles/ks 
of  products  of  detonation,  the  chemical  energy  Q released  in  detonation, 
and  the  actual  Cv  and  average  Cy  heat  capacities  can  be  calculated  for 
each  temperature  and  fugacity  factor  considered.  This  method  of  compu- 
tation is  given  in  appendix  II  of  T.  R.  V. 

The  method  of  solution  of  the  thermo-hydroclynamic  problem  is  similar 
to  the  method  given  in  the  second  part  of  T.  R.  VIII.  The  difference  is 
in  the  substitution  of  the  assumed  x ■ X(v2)  relationship  in  place  of 
the  empirical  velocity  vs.  density  relationship.  In  order  to  facilitate 
this  procedure,  the  following  expressions  have  been  computed  on  the  basis 
of  the  particular  selection  of  the  x **  X(^  mentioned  above. 

9o  ■ 1 * x - (f  2)  (4) 
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Table  I 

Corresponding  Relationships  of 
x,  v^,  and  log  F Under  the  Assumption  of  x ■ X(v2) 


X 

v2 

lo£,  v„ 

log10  F 

0 

o 

log10  F “ 

0.6 

2.050 

0.312 

0.35 

1.59 

0.36 

0.7 

1.736 

0.240 

0.57 

1.68 

0.58 

0.8 

1.514 

0.180 

0.78 

1.77 

0.79 

0.9 

1.343 

0.128 

1.01 

1.86 

1.02 

1.0 

1.202 

0.080 

1.24 

1.95 

1.26 

1.1 

1.085 

0.035 

1.48 

2.05 

1.50 

1.2 

0.992 

9.996- 

1.74 

2.15 

1.77 

1.3 

0.913 

9.961- 

2.02 

2.24 

2.05 

1.4 

0.845 

9.927- 

2.31 

2.32 

2.34 

1.5 

0.785 

9.895- 

2.62 

2.39 

2.66 

1.6 

0.732 

9.864- 

2.96 

2.45 

3.01 

1.7 

0.684 

9.835- 

3.33 

2.53 

3.39 

1.8 

0.643 

9.808- 

3.72 

2.61 

3.77 

1.9 

0.604 

9.781- 

4.15 

2.62 

4.25 

2.0 

0.568 

9.754- 

4.62 

2.62 

4.74 

2.1 

0.534 

9.727- 

5.13 

2.55 

5.30 

2.2 

0.499 

9.698- 

5.68 

2.51 

5.91 

2.3 

0.467 

9.669- 

6.29 

2.50 

6.61 

2.4 

0.437 

9.640- 

6.95 

2.50 

7.46 
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(5) 

(6) 

where 

N * log10  e * 0.4343 


log10F  = N<‘X  + J ~ “ log  v. 


K(v) 

NY  - N J'  (eX-l)|£ 
1 K(oo) 


and  Y^  is  part  of  the  integral  Y appearing  in  T.  R.  VI  and  VIII.  The 
values  of  these  expressions  corresponding  to  x*  s in  the  range  0.6,  0.7,..., 
2.4  are  given  in  table  1. 

To  calculate  the  thermo-hydrodynamic  quantities  for  various  conditions, 
the  equations  of  T.  R.  VIII  are  used,  and  the  following  method  is  suggested. 
Corresponding  to  a set  of  x values  over  a suitable  range,  say  0.6,  0.8,..., 
2.4,  values  for  T and  n are  selected.  Since  the  method  of  solution  is  one 
of  successive  approximations, the  choice  of  values  should  be  made  to  corre- 
spond to  that  which  is  expected.  The  integral 

NY  - N/n(v)(ex-l)^  (7) 

2 nW  n 

is  computed,  and  from  equation  (2.13)  of  T.  R.  VIII, 

log10  F = (lo^0  F - NYX)  + NY2  (8) 


Corresponding  to  the  sets  of  T and  log^g  F,  the  quantities  n,  Q,  C^, 
and  C are  obtained.  The  temperature  is  then  calculated  by 


where 


Q ♦ 


T,C 
1 Y 


C - nRe 
V 20 


(9) 


(10) 
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R - 1. 987*10  for  T in  degrees  Kelvin,  and  is  usually  taken  to  be 
300°K. 

If  the  calculated  and  assumed  values  of  T and  r.  agree,  then  these 
values  are  taken  as  part  of  the  solution.  In  general,  however,  a few 
iterations  will  be  required  to  satisfy  these  conditions.  It  is  suggested 
that  in  proceeding  from  one  iteration  to  the  next  the  calculated  values 
of  n be  used  from  the  previous  computation,  but  that  the  assumed  and 
calculated  values  of  T be  averaged.  This  method  was  found  to  be  success- 
fill  in  obtaining  rapid  convergence.  Also,  for  computing  the  integral  (7) 
and  the  derivative  in  (10), the  approximatio-  methods  given  in  T.  R.  VIII 
are  satisfactory. 

The  detonation  velocity  D and  initial  density  A^  are  calculated  by 
use  of  the  formulae 

D2  - nRT2eX  (1  ♦ 0)2/0  (ll) 


A 


1 


v2(i  ♦ e) 


(12) 


where  R * 8.314  for  velocity  in  meters/sec.  From  these  values  one  may 
provide  a plot  of  the  predicted  velocity  against  initial  density.  The 
other  thermodynamic  quantities  can  be  calculated  by  use  of  the  equations 
in  T.  R.  VI  or  VIII. 

As  an  illustration  of  this  method,  the  velocity  of  TNT  was  calculated. 
The  results,  designated  by  "Method  B",  are  compared  with  those  previously 
obtained  in  T.  R.  VIII,  "Method  A".  Table  II  gives  a tabulation  of  the 
thermo-hydrodynamic  variables,  and  figure  4 gives  a plot  of  the  velocity- 
density  relationship.  It  is  clear  that  the  results  of  Method  A and  B 
would  be  identical,  providing  the  x - X(v2)  relationship  had  been  selected 
as  that  calculated  by  Method  A.  The  difference  in  the  two  sets  of  values 
is  thus  the  measure  of  the  error  for  TNT,  introduced  by  the  fundamental 
assumption.  The  difference  is  insignificant  for  most  of  the  quantities. 

If  the  integrals  in  the  log  ? function  and  derivatives  in  the  0 
function  were  neglected  in  the  thermo-hydrodynamic  calculations,  it  would 
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be  possible  to  devise  a very  simple  procedure  for  predicting  detonation 
velocities.  Under  these  conditions  it  is  unnecessary  to  compute  a com- 
plete table  of  composition  against  temperature  and  fugacity  factors,  but 
sufficient  to  calculate  the  particular  composition  at  the  conditions  of 
interest.  Assuming  in  accordance  with  experience  that  the  velocity-den- 
sity relationship  is  linear,  solutions  of  the  thermo-hydrodynamic  problem 
at  two  specific  set  of  conditions  will  be  sufficient  to  determine  this 
relationship. 

The  following  method  was  found  satisfactory  in  solving  the  problem 
under  these  strong  assumptions.  For  a given  x,  obtain  from  table  I a 
corresponding  value  of  log  F.  Since  in  general  Y..  and  tend  to  cancel 
each  other,  to®  F ^ ? better  aDDrox'mation  for  the  fugacity  factor  than 
log  F - NY^.  Using  this  value  log  F,  and  assuming  a temperature  T2,  the 
corresponding  composition  and  n,  Q,  Cv,  and  C^,  are  calculated.  Substi- 
tuting these  values  into  equations  (9)  end  (10),  where  the  derivative 
in  (10)  is  neglected,  a value  for  T2  is  calculated.  Iterating  as  pre- 
viously suggested,  a consistent  solution  is  obtained  when  the  assumed  and 
calculated  values  of  agree.  The  velocity-density  relationship  is  then 
obtained  by  use  of  equations  (ll)  end  (12).  Two  such  solutions  determine 
the  desired  linear  relationship. 

As  an  illustration  of  this  method,  the  case  of  TNT  is  considered 
once  again.  The  results  are  given  in  table  II  under  Msthod  C,  and  in 
figure  4,  curve  C.  Once  again,  the  method  gives  results  which  agree  ap- 
proximately with  the  complete  and  accurate  solution.  Additional  calcu- 
lations verified  that  this  method  gives  a linear  velocity-density  rela- 
tion. Hence,  this  is  valid  assumption  for  this  scheme.  Furthermore,  the 
difference  between  curves  A and  B in  figure  4 indicates  the  significance 
of  a correct  functional  relation  x * X(v2).  The  difference  between  curves 
B and  C indicate  the  error  introduced  by  the  additional  assumptions  dis- 
cussed here.  Appreciating  the  magnitude  of  some  of  the  terms  neglected 
by  this  method,  one  can  only  explain  the  reasonable  results  on  the  basis 
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of  compensating  errors.  Incidentally,  for  computing  relative  velocities 
instead  of  absolute  ones,  the  approximation  Method  C is  actually  quite 
good. 

As  applications  of  this  method,  two  problems  were  considered  for  the 
case  of  TNT.  First,  the  effect  of  initial  temperature  or  detonation  velo- 
city, and  second,  the  effect  of  non-reaction  in  an  explosive.  Since  the 
calculations  are  made  by  this  approximation  method,  these  effects  are  to 
be  studied  qualitatively.  Corrections  could  be  made,  however,  by  extra- 
polation and  use  of  the  correct  velocity-density  relation. 

The  first  problem  is  answered  by  simply  varying  the  parameter  T^ 

(usually  taken  to  be  300°K)  in  equation  (9).  The  set  of  temperature 
values  200,  300,  LOO,  and  500°K  were  considered  for  the  case  of  TNT,  and 
the  velocity-density  curve  was  calculated  for  each  such  value.  The  results 
are  given  in  table  III.  The  effect  of  this  parameter  is  also  demonstrated 
in  figure  5,  where  the  velocity  density  curves  are  plotted.  It  is  to  be 
observed  that  the  effect  is  not  appreciable  on  ary  of  the  therj»-hydrodynamic 
quantities,  but  that  an  increase  of  initial  temperature  would  tend  to  in- 
crease the  detonation  velocity. 

The  second  problem  is  solved  by  considering  that  a specific  weight 
percentage  of  explosive  reacts.  That  is,  if  M^,...,  is  the  initial, 
composition  (in  gram/atom/lOOgram)  of  the  m element?  present  in  the  com- 
position, then  M^, . . . , N^M  is  the  composition  of  an  explosive  equiva- 
lent to  an  N reaction  of  the  first  explosive.  Calculations  were  carried 
w 

ouu  for  TNT  at  two  different  densities,  taking  the  set  0.60,  0.60,  and 
0.40  as  the  values  for  N^.  These  results  are  given  in  table  IVA,  and  the 
velocity-density  curves  are  given  in  figure  6.  The  effect  of  non-reaction 
is  clearly  demonstrated  by  these  results. 

In  T.  R.  IV  the  simple  approximation  equation 

-(2,)1/2 : (»M)1/2  (13) 
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Table  IV  A 


Effect  of  Fartial  Reaction  on  Relative  Velocity  of  TNT 


N 

w 

n 

Q 

Cv 

T*10~3 

A1 

d/d* 

1.00 

30.1 

945 

.340 

3.58 

1.04 

1.00 

0.80 

24.2 

755 

.284 

3.46 

1.04 

0.88 

0.60 

18.2 

544 

.204 

3.42 

1.04 

0.76 

0.40 

12.1 

350 

.132 

3.40 

1.04 

0.62 

1.00 

23.4 

1175 

.379 

4.02 

1.56 

1.00 

0.80 

18.7 

927 

.304 

3.99 

1.56 

0.89 

0.60 

13.9 

675 

.227 

3.95 

1.56 

0.76 

0.40 

9.1 

417 

.143 

3.93 

1.56 

0.61 

Table  IV  B 

Comparison  of 

N 1/2 

(n/n*)l/2. 

and  D/D 

for  TNT 

A1  ‘ : 

1.04 

A.  - 1.56 

N 

w 

Ni/2 

w 

(n/n*)^ 

D/D* 

<n/nV/2 

D/D* 

0.80 

0.894 

0.897 

0.88  3 

0.894 

0.889 

0.60 

0.775 

0.777 

0.760 

0.771 

0.762 

0.40 

0.632 

0.634 

0.624 

0.624 

0.613 
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FIG. 6 EFFECT  OF  NON-REACTION  ON  VELOCITY-DENS 
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was  derived,  where  D and  n are,  r3spectively,  the  detonation  velocity  and 
the  total  number  of  moles/kg  of  explosives  .products  for  the  incompletely 
reacting  explosive,  and  D and  n are  the  corresponding  quantities  for 
complete  reaction  under  'ideal'  conditions.  Table  IVB  gives  a tabulation 
of  the  three  terms  of  (13)  for  TNT  at  the  two  densities  and  set  of  values 
of  reported  in  table  IVA.  Since  ratios  only  are  here  considered,  the 
results  are  quantitative,  and  agree  with  the  prediction  of  equation  (13). 

Detonation  Velocity  in  Explosives  with  Solid  Inert  Additives 

As  is  clear  from  the  study  of  the  preceeding  section,  air  is  an  effect- 
ive diluent  for  explosives.  The  velocity-density  relationship  i3  practi- 
cally linear,  with  velocity  decreasing  as  the  volume  fraction  of  air  is 
increased.  It  is  of  interest  to  study  the  effect  of  other  diluents  in 
pure  explosives.  The  following  discussion  is  concerned  with  the  effect  of 
adding  solid,  non-explosive  materials  to  an  explosive. 

Let  M and  A^  be  the  mass  and  density,  respectively,  of  the  explosive- 
inert  mixture,  and  M.,  and  A,  the  correspondir^  quantities  of  the  component 

i,  (pure  explosive  or  solid  inerts).  The  crystal  density  A of  the  mixture 

c 

is 

Ac  - (I  M.)(Z  ^r1  - Ax  (1  - A)”1  (14) 

where  i runs  over  the  indices  of  the  solid  componerts,  and  A is  the  volume 
fraction  of  air  in  the  mixture.  That  is, 

Ax  - Ac  (1  - A)  (15) 


ai 

{ 


1 


u 

i 

i 

l 

l 

•* 

A 


In  the  following,  the  same  volume  fraction  of  air  will  be  assumed  in  the 
mixtures  as  in  the  case  of  the  ideal  pure  explosive.  That  is 


A 


1 


(16) 


I 
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where  A^  is  the  crystal  density  of  the  pure  explosive.  Here  only  one 
solid  diluent  is  considered.  Hence,  under  these  conditions,  if  N^  is  the 
weight  fraction  of  explosive,  1 - N^  is  the  weight  fraction  of  the  inert, 

and  N 1 - N -1 

4c  - ^ ‘-17s  > <17) 


where  A_  and  AT  are  the  crystal  densities  of  the  explosive  and  inert, 
respectively. 

In  addition  to  these  preliminary  considerations,  the  following  were 
applied  in  calculation  of  the  detonation  velocity  of  an  explosive-inert 

mixture.  Using  the  covolume  equation  of  state  (equation  (2))  the  ideal 

* 

detonation  velocity  D is  given  by 


* * 
Ti-° 


V4  ("*^*)1/2 

(9  )1/2  2 


where 


* * 
* Cy  + n R 


* 

- (— ) 

Ws 


Here  the  relationships  (18)  and  (19)  refer  to  the  ease  of  a pure  ideal 
explosive,  and  corresponding  non- starred  equations  describe  the  ideal 
conditions  of  an  explosive-solid  inert  mixture. 

For  all  practical  purposes  if  the  explosives  are  compared  at  the  3ame 

v / * \ //  #\l/2 

relative  free  space,  8 ■ 3 . Furthermore,  the  term  O ♦ l)/(p  ) is 

* 

quite  insensitive  to  0 for  the  usual  range  of  values  of  this  quantity. 

* * 

For  example,  a 20  per  cent  increase  in  the  value  of  0 near  0 “0.5 

corresponds  to  less  than  3 per  cent  decrease  in  the  value  of  (0  ♦ 1)/ 

(0  ) ' . Hence,  by  dividing  the  unstarred  relationship  corresponding  to 
(18)  by  the  equation  (18), 

r a nTy  1/2 

K - (-H;)  ‘ (20) 

D a n 
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where 

a - (1  - A^)"1,  a*  - (1  - A*  a*)"1  (21) 

The  situation  of  an  explosive-inert  mixture  is  similar  to  the  problem 

of  non-reaction  in  a pure  explosive.  Here,  unreacted  explosive  is  to  be 

considered  as  inert.  Under  these  conditions,  as  is  indicated  on  page  20 

* 

of  T.  R.  VI,  T£  * T2.  This  is  demonstrated  also  by  the  results  reported 

in  table  IVA  for  the  case  of  TNT.  At  A^  ■ 1.56,  the  error  is  less  than 

2.5  per  cent  in  assuming  the  complete  reaction  temperature  for  the  case 

of  40  per  cent  reaction  in  pur*  TNT.  Furthermore,  assuming  complete  re- 

♦ 

action  of  the  pure  explosive,  n/n  ■ N^.  Therefore,  from  (20), 

”*  “ "*  (»J1/2  (22) 

Da 

If  the  solid  inert  is  considered  to  be  non-reacting  explosive,  then 

it  is  clear  by  (21)  that  equation  (22)  would  reduce  to  (13).  For  under 

these  conditions,  the  initial  density  and  co-volume  of  the  ’’mixture"  is 

the  same  as  the  pure  explosive. 

* * 

Assuming  that  D and  a are  known  for  a particular  ideal  explosive, 

D can  be  calculated  for  a specific  in  terms  of  a.  Also,  by  (21),  a 
is  a simple  function  of  the  density  A^  and  co-volume  a of  the  mixture. 

Since  A^  is  readily  obtainable,  the  only  problem  is  in  obtaining  a. 

To  calculate  the  co-volume  of  the  mixture,  the  following  is  considered. 
First, 

a - Og  + Oj  (23) 

where  ou  is  the  part  of  the  total  co-volume  contributed  by  the  explosive, 
and  cij  that  of  the  inert.  It  is  assumed  that  * N^a  (p)  where  a (p) 
is  the  co-volume  of  the  ideal  pure  explosive,  corresponding  to  the  pres- 
sure p of  the  mixture  explosive.  Furthermore,  ■ (1  - N^)  v^  (p),  where 
vT  (p)  is  the  specific  volume  of  the  inert  at  pressure  p. 

* 4. 

The  parameter  a (p)  is  obtained  from  the  solution  of  the  thermo-hydro- 
dynamic  problem  of  the  pure  explosive.  It  is  simply  the  value  of  a 
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corresponding  to  the  pressure  p.  If  the  complete  co-volume-pre3sure  re- 
lationship has  not  been  calculated  for  a particular  explosive,  the  follow- 
ing approximation  is  satisf -itory.  From  hydrodynamic  theory 

(24) 


^DW 


and 


D 

W 


V — V 
1 2 


(25) 


If  it  is  assumed  now  that  v ^ ■ cv^,  where  c is  constant,  then 

p - (1  - c)D2 


(26) 


If  Vg  is  known  for  a particular  A^,  then  c = v^A^.  This  assumption  is  a 
valid  assumption  for  most  solid  explosives,  and  corresponding  values  of 

v„  and  A,  are  obtainable  for  various  important  explosives  in  "Theory  of 

2 1 (1) 

Detonation"''  . 

As  part  of  the  procedure  for  calculating  the  velocities  of  the  explo- 

. . 

sive-inert  mixtures,  it  is  advisable  that  a plot  be  made  of  a (p;  vs. 

. * * 

p/p  , where  p is  the  pressure  of  the  pure  explosive,  corresponding  to  a 
* 

density  A^.  3uch  a plot  is  given  in  figure  7 for  TNT  at  two  different 
densities.  Similar  curves,  which  parallel  these,  are  easily  obtainable 
for  other  explosives. 

The  calculation  of  v^  (p)  depends  upon  the  following  properties  of 

solids.  It  is  convenient  to  use  the  identity 

vT(p) 


i (p)  =J. 


(27) 


where  v^  is  the  specific  volume  of  the  inert  at  one  atmosphere  pressure, 
and  Aj  is  the  initial  density  of  the  inert . Introducing  the  compressibi- 
lity factor  3 (p)  of  the  inert  at  pressure  p. 


Vj(p) 


■ 1 - 30(p)p 


1 
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Also  in  "Properties  of  Solids 


„(2) 


Cook  has  derived  the  following  rela- 


tionship between  the  compressibility  at  atmospheric  pressure  and  that  at 


pressure  p 


1 - eQp  ♦ o0p r 

0o(p)  “ 0o  (1  ♦ 2/3  3°p)(l  * 1/3  0QP) 


(29) 


where  3 

o 

since 


is  the  compressibility  of  the  solid  inert  at  one  atmosphere.  Also 


* 

P 


* * 
D 


nRT2  1/2 


(30) 


under  the  previous  assumptions, 


£ 

* 

P 


-*  (V1/2 

D 


(31) 


The  procedure  in  calculating  D for  specific  values  of  N , A.  and  3 , 

^ ^ W JL  U 

and  known  D , p , and  A.,,  is  as  follows.  Assuming  a value  of  the  ratio 
D/D  , the  ratio  p/p  is  computed  by  means  of  equation  (31).  The  value  of 
3Qp  is  then  obtained,  and  by  (29)  the  value  of  3q(p)p,  hence,  by  (28)  the 
value  of  vt(p)/yt.  The  value  of  the  co-volume  a is  iiunediately  obtained, 

using  the  density  of  the  inert,  the  a (p)  vs.  p/p  relationship,  and 

# 

equation  (23).  Hence,  by  (21)  and  (20),  a calculated  value  of  D/D  is 
obtained.  If  this  computed  value  agrees  with  the  assumed  ratio,  the  so- 
lution is  complete.  If  it  does  not  agree,  then  the  process  is  repeated, 
using  the  average  value  of  the  assumed  and  calculated  values  of  D/D  . This 

iteration  procedure  usually  converges  very  rapidly. 

(3) 

NOL  reported  a set  of  experimental  data  for  TNT,  mixed  respectively 

* 

with  NaCl,  Ba(N0^)2,  PbSO^,  and  Cu.  A value  of  D * 6670  raeters/sec.  was 
reported  here  for  a density  « 1.59.  This  'ideal* velocity  value  wag  used 
in  thi s study,  even  though  it  is  low,  since  the  calculated  and  observed 
velocities  are  compared  directly.  Also,  for  this  study  the  experimental 
velocities  are  corrected  for  the  same  relative  air  space  in  the  mixtures 
as  in  the  pure  TNT  measurements.  The  results  of  this  calculation  are 
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shown  in  table  V.  Figure  3 is  included  in  the  case  of  TNT  - Cu  comparing 
graphically  the  experimental  and  calculated  results. 

Further  calculations  were  made  for  TNT  - NaCl,  TNT  - Al,  TNT  - NaNO^> 
50-50  Per.tollte  - NaCl,  RDX  - NaCl,  and  Comp  B - Al,  assuming  conditions 
for  no  reaction  of  the  Al  and  NaNOy  These  results  are  given  in  table  VI. 

The  velocity  data  for  the  pure  explosives  listed  here  were  taken  from  meas- 
urements of  NDRC  at  Bruceton^ . References  (l)  and  (5)  were  used  to  ob- 

* 

tain  the  compressibilities  of  the  inert  substances.  The  a vs.  p relation 
for  TNT  is  given  in  figure  7j  the  same  relation  for  the  other  three  explo- 
sives was  obtained  by  the  above  mentioned  approximation  method,  using 
values  for  c calculated  from  the  data  given  in  "Theory  of  Detonation"^ 

Figures  9 and  10  are  included  to  indicate  the  effect  of  dilution  on  detona- 
tion velocity  over  a range  of  densities. 

Experimental  work  is  now  in  progress  on  this  project  to  obtain  velocities 
for  some  of  the  mixtures  considered  above.  Wherever  available,  these  re- 
sults are  reported  in  table  VI.  A comparison  of  the  calculated  and  experi- 
mental values  reveals  a satisfactory  agreement. 

It  is  to  be  noted  that  if  a measured  velocity  of  a particular  mixture 
is  given,  then  by  an  'inverse'  method  of  solution,  the  compressibility  of 
the  inert  can  be  calculated.  This  would  permit  one  to  'measure'  compresi- 
bility  of  solids  at  pressures  considerably  beyond  those  reported  by  Bridgman. 
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(3\ 

Table  V:  Comparison  of  Calculated  and  Observed 

Velocities  for  TNT-Inert  Mixtures 


TNT; 

D*  - 6670* 

- 1*59 

• 

Calculated 

Observed* 

Inert 

■w 

*1 

D-IO-3 

D-10"3 

Nad 

0.90 

1.61 

6.52 

6.54 

0.65 

1.71 

6.11 

6.20 

0.50 

1.78 

5.85 

5.67 

Ba(N03)2 

0.90 

1.65 

6.41 

6.53 

0.50 

2.08 

5.30 

5.48 

0.30 

2.39 

4.71 

4.66 

0.27 

2.44 

4.62 

4*60 

Cu 

0.90 

1.73 

6.47 

6.63 

0.80 

1.90 

6.22 

6.31 

0.70 

2.10 

5.94 

6.12 

0.60 

2.36 

> 5.66 

5.56 

0.50 

2.68 

5.37 

5.26 

0.40 

3.12 

5.07 

5.25 

0.30 

3.70 

4.76 

4.86 

0.27 

3.92 

4.67 

4.70 

PbSO, 

4 

0.90 

1.72 

6.45 

6.53 

0.80 

1.86 

6.21 

6.35 

0.70 

2.04 

5.96 

6.15 

0.50 

2.52 

5.42 

5.60 

0.30 

3.27 

4.79 

4.87 

O.Z' 

3.42 

4*68 

4.8S 

^ NavOrd  Report  1538.  Velocities  corrected  to  same  air  space  as  pore 
TNT  measurements.  D*  * 6670  is  low  for  the  ideal  Telocity  at  this 
density.  It  was  used  uncorrected  since  only  the  ratio  of  D/D*  is 
of  interest. 
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Table  71;  Velocity  Predictions  and  some  Observations 


Mixture 

* 

A. 

D*  *10~3 

N 

A. 

D*10-3 

D-10-3 

1 

V 

1 

(Calc) 

(Obs.) 

TNT-NaCl 

1.59 

6.91 

0.90 

1.63 

6.77 

6.75 

1.59 

6.91 

0.70 

1.71 

6.49 

6.51 

1.59 

6.91 

0*60 

1.76 

6.34 

6.27 

TNT-A1* 

i 

1.56 

6.81 

0.80 

1.69 

6.52 

TNT-NaNO^* 

1.59 

6.91 

0.50 

1.84 

6.18 

RDX-NaCl 

1.26 

6.86 

0.70 

1.32 

5.96 

5.80 

1.36 

7.20 

0.40 

1.51 

4.72 

5.00 

Pentolite-NaCl 

1.60 

7.34 

0.70 

1.83 

6.94 

6.93 

1.20 

6.10 

0.70 

1.28 

5.31 

Comp  B-Al* 

1.65 

7.64 

0.80 

1.78 

7.55 

* Assuming  no  reaction  in  the  detonation  head.  This  should  apply 
approximately  foi  65  ■»  mesh  A1  and  20-30  mesh  fiaNO.  for  charges 
in  5 cm.  diameter. 
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